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^^ ' Abstract 

f~^ I In the presence of a static quark-antiquark pair, the spectrum of the low-lying states in 

SU(A'") gauge theories is discrete and likely to be described, at large quark separations r, 

by an effective string theory. The expansion of the excitation energies in powers of 1/r, 

O ,' which derives from the latter, involves an increasing number of unknown couplings that 

4j ' characterize the string self-interactions. Using open-closed string duality, we show that the 

possible values of the couplings are constrained by a set of algebraic relations. In particular, 

the corrections of order 1/r^ must vanish, while the 1/r^ terms (which we work out for the 

^^ . few lowest levels) depend on a single adjustable coupling only. 
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1. Introduction 

The formation of string-like flux tubes in SU(A^) gauge theories is not a proved fact, 
but the results obtained in lattice gauge theory over the years leave little doubt that 
this physical picture is basically correct. Recent high-precision numerical simulations 
now even confirm the presence of characteristic string fluctuation effects in the static 
quark potential [1,2] and in ratios of loop correlation functions [3-7]. 

A puzzling aspect of these findings is, however, that string behaviour appears to 
set in at quark-antiquark separations r as small as 0.5 to 1 fm or so, where the energy 



distribution in the ground state of the gauge field is stih far from being tube-hke. 
Moreover, at these distances the spectrum of the low-lying excited states is poorly 
matched by the regular pattern (with level splittings equal to vr/r) that is obtained 
in the case of a free string with fixed ends [8,9,2]. Closed flux tubes wrapping around 
a compact space dimension have also been considered, with similar conclusions [10- 
14]. 

At present none of these results exclude that the observed mismatches will grad- 
ually disappear at larger distances r. The energy spectrum should in any case be 
compared with the predictions of the effective string theory of refs. [15-17], in which 
the string self-interactions are taken into account in a systematic manner. According 
to this theory, the free-string energies are only the leading terms in an asymptotic 
expansion in powers of 1/r, while the higher-order corrections arise from interaction 
terms of increasing dimensions. There are finitely many terms at each order and a 
corresponding number of effective coupling constants. 

Our principal goal in this paper is to work out the expansion of the few lowest 
energy levels up to order 1/r^. We shall first show, however, that the a priori un- 
known values of the couplings in the effective theory must satisfy certain algebraic 
relations in order to be consistent with the closed-string interpretation of the parti- 
tion function. As a consequence, corrections proportional to 1/r^ are excluded and 
the number of independent couplings at the next order of the large-r expansion is 
reduced from 2 to 1. The predictive power of the effective theory is thus significantly 
enhanced. 



2. Eflfective string theory 

In this section we define the effective string theory and explain which properties of 
the gauge theory it is expected to describe. We also compute the string partition 
function to next-to-leading order, partly for illustration and partly because this will 
serve as a starting point for the discussion in the following section. Much of what 
is being said here already appeared in ref. [1], and we refer the reader to this paper 
for a less abbreviated presentation. 

2.1 Definition 

Eventually we are interested in the string energy spectrum in 4 space-time dimen- 
sions, but string formation is also studied in 3 dimensions and we therefore consider 



the effective tfieory in any dimension d > 3. Tlie ffuctuations of a string witli fixed 
ends can be parametrized by a displacement vector h{z), with d — 2 components in 
the directions orthogonal to the string axis, where 

z={zo,zi), 0<zi<r, (2.1) 

are the world-sheet coordinates. Fixed ends means that the field satisfies 

M^)L,=o=M^)L,=. = (2.2) 

for all times zq, and in the time direction we shall impose periodic boundary condi- 
tions, with period T, at least as long as we are concerned with the string partition 
function. 

The action of the effective theory is written in the form of a series 

S = arT + fiT + So + Si + ---, So = ^ f d^z (dahdah) , (2.3) 

in which a denotes the string tension, ^ an arbitrary mass parameter and Si, S2, ■ ■ ■ 
string self-interaction terms of increasing dimension. For the calculation of the string 
energy levels it is convenient to choose the metrics in space-time and on the world- 
sheet to be euclidean. We shall use dimensional regularization and thus extend the 
world-sheet to D = 2 — e dimensions at this point (see appendix A). 

As is generally the case in field theories on manifolds with boundaries [18,19], the 
possible interaction terms are localized either on the world-sheet or at its boundaries. 
They may only depend on the derivatives of h{z) and must respect the obvious 
symmetries of the system. Moreover, any terms that are formally removable by 
a field transformation can be dropped, since their effects on the energy spectrum 
amount to a renormalization of the coupling constants multiplying the interactions 
of lower dimension. 

Inspection now shows that 

Si = \h j d^-'z{{dihdih)^^^^ + {dihdih)^^^^] (2.4) 

is the only possible interaction term with coupling b of dimension [length] . At second 
order the couplings have dimension [length]^ and the complete list of terms is 

^2 = \c2M' I d'^z (dahdah) (dbhObh) , (2.5) 
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\c^M' I d^z {dahdbh) (dahdbh) , (2.6) 



where M denotes the renormahzation scale and the brackets indicate in which way 
the space indices are contracted. Since we are using dimensional regularization, there 
are no power divergences and b is thus not renormalized. In principle an additive 
renormahzation of C2,C3 proportional to b^ may be required, but it turns out that 
there are no logarithmic divergences to this order. 

2.2 Partition function and the Polyakov loop correlation function 

In the effective theory the partition function and the correlation functions of the 
field h{z) can be expanded in powers of the couplings 6, C2, C3, . . . following the usual 
steps. From this point of view the theory looks just like any other non-renormalizable 
quantum field theory. In particular, the perturbation series really is a low-energy 
expansion, where the higher-order terms are proportional to inverse powers of the 
external distances such as T and r. 

To leading order the exact expression for the partition function 

Zo = e-"'-^-^^r/(g)2-'', (2.7) 

V(.q)=q^'ll(^-qn, <? = e-"^/^ (2.8) 

ra=l 

is known since a long time (see refs. [20,21,3] for example; the index "0" reminds 
us that we are dealing with the free-string theory). When expanded in powers of q, 
this formula leads to the representation 

00 
Zo = Y,Wne-''"^, (2.9) 

n=0 

with some positive integral weights Wn and energies E^ given by 

S0=ar + /i + -{-^(d-2)+n}. (2.10) 



In particular, the splitting between successive energy levels is equal to vr/r to this 
order of the perturbation expansion. 

In the presence of a static quark-antiquark pair, the spectrum of the low-lying 
states in the (pure) SU(A^) gauge theory is surely more complicated than this, but 
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it is conceivable that the effective string theory provides an asymptotically correct 
description of the energy levels at large distances r. The partition function in this 
sector of the gauge theory is proportional to the Polyakov loop correlation function 
{P{x)*P{y)) at distance r = |x— y| on a euclidean space-time cylinder with circum- 
ference T in the time direction [1]. Evidently, if the energy spectra match (including 
multiplicities), the same must be true for the associated partition functions. 

We would like to emphasize, however, that the (full) string partition function Z 
is only expected to match the Polyakov loop correlation function to any finite order 
of the perturbation expansion in the effective theory, i.e. only at large T and r and 
up to higher-order terms. Decays of highly excited states through glueball radiation 
are in fact not included in the string theory, at least not obviously so, and this must 
be reflected by relative differences in the partition functions of order e"™""^, where 
m denotes the mass of the lightest glueball. 

2.3 Partition function at next-to-leading order 

The order at which a given interaction term must first be included in the perturba- 
tion expansion coincides with the dimension of the associated coupling constant. At 
next-to-leading order the partition function thus reads 

Z = Zo{l-(5i)o}, (2.11) 

where the expectation value (. . .)o is to be evaluated in the free-string theory. If we 
introduce the free propagator, 

{hkiz)hiiw))Q = 6kiGiz,w), (2.12) 

the first-order correction term assumes the form 

(5i)o = \h{d-2) j d^-^z { aiG(z, z) ai i.^^o + dMz, ^)^i L^=,}, (2.13) 

and the rules of dimensional regularization must now be invoked to work it out. 

Deferring the details of the calculation to appendix A, we note that there are no 
ultra-violet divergences and quote the result 

(5i)o = -6(d-2)|-ln,7((7). (2.14) 

The correction thus amounts to a shift r ^ r — 6 in the second factor of the leading- 
order expression (2.7) (this simple rule actually extends to the next order in 6, but 
we shall not need to know this here). 



3. Open closed string duality and its consequences 

We now proceed to show that the values of the coupling constants 6, C2, C3, . . . in the 
effective string theory are constrained by certain algebraic relations if the Polyakov 
loop correlation function is to be matched by the string partition function at large 
T and r. 

3.1 Dual interpretation of the Polyakov loop correlation function 

Usually the Polyakov lines are considered to run along the time axis, forming closed 
loops through the periodic boundary conditions in this direction. This leads to the 
interpretation of the Polyakov loop correlation function as the partition function of 
the gauge theory in the presence of a static quark-antiquark pair. 

We may, however, just as well take the compact dimension to be space-like, in 
which case the field-theoretical interpretation of the Polyakov loop P{x) is that it 
creates (or annihilates) a closed fiux tube wrapping around the world at, say, time xi 
and transverse spatial position x± = (x2, • • • .,Xd-i)- The loop correlation function 
at zero transverse momentum is therefore expected to have a spectral representation 
of the form [22-25] 

/oo 
dx^(P(x)*P(y)) = j;|t;„|%-^"l--^^l, (3.1) 

n=0 

where the exponents En and the coefficients Vn are the energies and transition matrix 
elements of the possible intermediate states f. As shown in appendix B, this implies 
that the correlation function itself can be expanded in a series of Bessel functions 
according to 

{PixyPiy)) = Y^\v^\^2r[^\ K.(,_3)(^„r). (3.2) 

n=0 \ / 

While the energy values and the transition matrix elements are not known, this equa- 
tion (and similarly the spectral representation in the open-string channel) severely 
constrain the possible analytic form of the Polyakov loop correlation function. 



f As in the case of the flux tubes with flxed ends, the higher excited states are unstable against 
glueball radiation and the spectrum of these states is thus continuous. This subtlety can be ignored 
in the following, since we shall be exclusively interested in the low energy values. 



3.2 Open-closed string duality at leading order 

An interesting question is now whether the partition function in the effective string 
theory can be expanded, at large T and r, in a series of Bessel functions of the form 
(3.2). If not this would imply that the effective theory cannot possibly provide an 
asymptotically exact description of the Polyakov loop correlation function. 
Using the modular transformation property of the ?7-function, 

r/(g) = (^) %(g), q = e-^-^/^, (3.3) 

the leading-order expression (2.7) for the partition function can be written as a series 
of exponentials, 

Zo = e-^^(^] l:u;„e-^"^ (3.4) 
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n=0 



where the weights Wn are the same as in the expansion (2.9) of the partition function 
in the open-string channel. Evidently, the closed-string energies 

zLtt" 
^°=aT+-{-^(d-2)+n} (3.5) 

must coincide with the energies En in eq. (3.2), up to terms of order b/T'^, if the 
series (3.2) and (3.4) are to match at large T and r. 

Since their arguments grow proportionally to aTr in this limit, the Bessel func- 
tions in the series (3.2) may be expanded, 



(PixyPiy)) = y:Kr{i^J e— { 1 + ^'^ ^l^'^^ ^^ +•••}, (3.6) 

where, to leading order, we may replace En by En and drop all subleading terms in 
the curly bracket. The free-string partition function (3.4) is then seen to be of the 
required form, in any dimension d, and open-closed string duality thus holds to this 
order of the perturbation expansion. 




3.3 Duality at next-to-leading order 

The only interaction term that contributes at next-to-leading order is the boundary 
term 5*1, and the partition function to this order is given by 

Z = (1 - aTh) Zo - 6^Zo 



^) j:-„|l + 6(^0-ar) + A(d_2)|e-^> (3.7) 



(cf. subsection 2.3). From the discussion of the leading-order expression of the par- 
tition function we know that E^ = E^ + 0(6), and once this relation is inserted in 
eq. (3.6) it becomes clear that the b/r term in the curly bracket in eq. (3.7) cannot 
be matched. We thus conclude that 

6 = (3.8) 

is the only possible value of this coupling, which is consistent with the requirement 
of open-closed string duality. 

3.4 Duality at higher orders 

At the next order of perturbation theory there are two string interaction terms, 5*2 
and 5*3, with couplings C2 and C3. Their contribution to the partition function is 
discussed in section 4, where we shall find that open-closed string duality holds if 
and only if 

{d-2)c2 + C3 = ——. (3.9) 

2a 

The number of adjustable couplings is thus reduced to 1, and actually to in di- 
mension d = 3, since the interactions S2 and 5*3 are proportional to each other in 
this case (cf. section 6.2). 

It is quite clear that the requirement of duality will lead to similar constraints at 
any order of perturbation theory. While this cannot be rigorously excluded, it seems 
unlikely that too many equations will be obtained at some point, because there are 
exact string partition functions with the correct properties in both the open and 
the closed string channels. The formal canonical quantization of the Nambu-Goto 
string, for example, yields such a partition function (appendix C). 



4. Partition function at t^vo-loop order 

We now always set 6 = so that the first non-trivial corrections to the free-string 
partition function come from the quartic interactions 5*2 and S3. At two-loop order 
the partition function is then given by 

Z = Zo{l- {82)0 -{83)0}, (4.1) 

where, as before, Zq and the brackets (• • •)o stand for the partition function and the 
expectation values in the free theory. 

4.1 Computation of {82)0 and (6*3)0 

The contraction of the fields in the interaction terms S2 and S3 results in tadpole 
diagrams with two loops attached to the interaction vertex. We evaluated them using 
dimensional regularization, as described in appendix A, but the same diagrams were 
actually already calculated many years ago by Dietz and Filk [20] who employed a 
(^-function finite-part prescription. 
In terms of the series 

Si(g) = q- lnr?(g) = _ - ^ ng" (1 - g")"^ , (4.2) 

n=l 



dq^ 



n=l 



(which are closely related to the first and second Eisenstein series) , the result of the 
computation reads 

{82)0 = C2 (d- 2) ^ {id-l)52iq) + 25,iqr} , (4.4) 

{83)0 = C3 (d - 2) ^ { 1 (d + 1) 52{q) + {d-l) Si(g)2} . (4.5) 

In particular, the partition function is again a power series in q = e~'^'^'^ , as in the 
case of the free-string theory. 



J^.2 Energy spectrum 

Once the string interactions are switched on, the accidental degeneracies of the en- 
ergy eigenstates, which are characteristic of the free-string theory, are not guaranteed 
to persist. The interpretation of the partition function as a spectral sum must take 
this possibility into account, and we thus write 

OO kn 

^ = EE^-.^«"''""^' (4.6) 

n=0 i=l 

where kn is the number of different energy values En^i that are equal to E^ to lowest 
order. Evidently the associated multiplicities Wn^i must satisfy the sum rule 

^Wn,i=Wn (4.7) 

i=l 

for all n > 0. 

It follows from this discussion that the partition function at two-loop order is of 
the general form 

OO 

Z = Y,w^{l-[Ei]T}e-^'^^, (4.8) 

n=0 

where [E^] is a linear combination of C2 and C3. This is just the weighted average, 
with weights equal to Wn,i/wn, of the first-order energy shifts E^ ^ in the multiplet 
with principal quantum number n. In particular, 

[El] = [El] + ^ {(12d - 14) C2 + (5d + 7) C3} , (4.10) 

and the average energy shifts of the higher levels can be obtained straightforwardly 
by expanding the partition function in powers of q. Note that the degeneracies of the 
ground state and of the first excited states are unchanged by the interactions. There 
is thus a single level in these two cases, with first-order energy shift -E^ ^ = [El^. 

Using the fact that the series Si and S2 are derivatives of the free-string partition 
function, the result 

[Ei] = [El] + ^n[^{d-2)-n]c2 if C3 = -2c2, (4.11) 

10 



is now also easily established. Couplings related in this way are a rather special case 
in many respects. In particular, the techniques introduced in section 5 can be used 
to show that the free-string degeneracies of the energy levels are preserved by the 
interactions at this order of perturbation theory. 

4-3 Closed-string representation 

For simplicity we set /x = in the following, since the associated exponential factor 
only complicates the formulae, but is of no importance for the open-closed string 
duality of the partition function. The duality property of the free-string partition 
function may then be compactly expressed through 



Zc 



j.^i(d-2) 

2r 



Z, 



Olr^2r 



(4.12) 



where the notation on the right implies a substitution of T by 2r and simultaneously 
of r by r/2. 

Similar identities can be derived for the first-order perturbations (4.4) and (4.5) 
since these are obtained from the ry-function by differentiation. Explicitly one finds 
that 



{S2 + 5*3)0 — ('S'2 + S3)o\j,^2r 



+ [{d - 2) C2 + C3] (d - 2; 



Att 



Sl(^) 



T 



IGvrr 



(4.13) 



and the closed-string representation of the partition function thus reads 



z = 


(-) {"'- 


^2r 








+ [{d - 2) C2 + C3] 


"d-2 
4Tr 


- a - 


1 d 
Tdr 



z, 



0It^2t 



When expanded in powers of q, this formula yields the series 



W-^) 



n=0 ^ 

+ [{d - 2) C2 + C3] 



-Elr 



^-[El]r 



47r r 1 , , , . d-2 
-{-±{d-2)+n] + ^^ 



(4.14) 



(4.15) 
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where the average closed-string energy shifts are given by 

[El] =2 [El] ^^^f^. (4.16) 

As in the case of the open-string channel, it is not possible, however, to deduce the 
degeneracy pattern of the energy levels from these equations. 

The comparison with the right-hand side of eq. (3.6) now shows that open-closed 
string duality holds to this order of perturbation theory if and only if the couplings 
C2 and C3 are related to each other as in eq. (3.9). In d = 4 dimensions this implies, 
incidentally, that the couplings lie on the special line C3 = — 2c2, where the accidental 
degeneracies of the energy levels in the open-string channel are not lifted by the 
interactions. 



5. Computation of the energy shifts El j 

The first-order shifts of the individual energy levels can be extracted from the con- 
nected correlation functions of the field h{z) at large time separations (we set T = 00 
in this section). Exactly which correlation functions are considered should not mat- 
ter in principle, but the computations can be simplified by introducing fields that 
create and annihilate definite free-string modes. We now first explain this and then 
proceed with a discussion of the symmetry classification of the low-lying energy 
eigenstates and of the associated eigenvalue equation. 

5.1 Creation and annihilation operators 
Let us consider the fields 



afc(^o,p) 



/ dzi / Az2...dzD-isni{pzi){p±do)hk{z), (5.1) 

Jo Jo 



where p = nir/r, n G N, is an arbitrary spatial momentum and k = 1, . . . ,d — 2 the 
space index that labels the directions transverse to the string axis. Since we intend 
to use dimensional regularization, the world-sheet is extended to D dimensions, as 
before, and an integration over the extra coordinates is included in eq. (5.1). It is 
straightforward to show that 

{at{zo,p)hi{w))^ = 6ki9{wo - zo)e-P^''"-'"^ sm{pwi), (5.2) 

{a^{zo,p)hi{w))^ = 6kie{zo - wo)e-P<^'"-^-^ sm{pw,), (5.3) 
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which suggests that a'^{zQ,p) and a^(zo,p) create and annihilate harmonic string 
modes at time zq with wave number n = pr/n. The normahzation of the associated 
one-mode states can be inferred from 

{a^{zo,p)a+{wo,q))^ = 5u6p^e{z^ - w^)e-P^'^-^^\pL-\ (5.4) 

and we also note that 

(a+(2;o,p)a+(wo,9)>o = {H{^o,p)ai {wo,q))^^ = 0, (5.5) 

which is entirely in line with the interpretation oi a'^{zQ,p) and a'^{zo,p) as creation 
and annihilation operators (for simplicity equal-time contact terms were dropped in 
these equations). 

Products of these fields create and annihilate states with any number of excited 
string modes. Effectively we have the association 

al^{zo,pi)...al^{zo,Pm) < — > |ni, . . . ,n^; fei, . . . , A;^) (5.6) 

between products of creation operators and the states in the free-string Fock space 
with transverse space-time indices fci , . . . , km and wave numbers 

nj = ^-, j = l,...,m. (5.7) 

Wick's rule ensures, incidentally, that the scalar products of these states are correctly 
given by the correlation functions of the corresponding products of the a"*" and a~ 
fields (the fields should be inserted at times infinitesimally different from each other 
so that contact terms do not arise). 

5.2 Symmetry considerations 

The Fock states (5.6) transform as tensors of rank m under the group 0{d — 2) of 
rotations about the string axis. Since the effective string theory is rotation- invariant, 
the energy eigenstates can be classified according to the irreducible representations 
of this group. 

For general dimension d and principal quantum number n = ni + - • ■ + njn < 3, the 
Fock space decomposes into altogether 11 multiplets, as shown in table 1. The states 
listed in the table are Fock states, or linear combinations of these, that transform 
according to some irreducible representation of 0((i — 2). As usual the ground state 
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Table 1. Energy eigenstates with quantum number n < 3 



n,i 


State 


0(d - 2) multiplet 


Wn,i 


0,1 


|0) 


scalar 


1 


1,1 


\l;k) 


vector 


d-2 


2,1 


|1,1) 


scalar 


1 


2,2 


\2;k) 


vector 


d-2 


2,3 


\l,h{k,l}) 


symmetric 2-tensor 


\d{d - 3) 


3,1 


|1,2) 


scalar 


1 


3,2 


|3;A;) 


vector 


d-2 


3,3 


|l,l,l;fe) 


vector 


d-2 


3,4 


\l,2;{k,l}) 


symmetric 2-tensor 


\d{d - 3) 


3,5 


\l,2;[k,l]) 


antisymmetric 2-tensor 


i(d-2)(d-3) 


3,6 


\l,l,l;{k,l,j}) 


symmetric 3-tensor 


i(d-2)(d + 2)(d-3) 



is denoted by 1 0) and the states that are non-trivial combinations of Fock states are 
given by 



d-2 



1,1)= J^|l,l;fc,fc), 



fc=i 



'kl 



d-2 



11,1) 



l,l;{A:,O) = |l,l;fc,0 



l,2) = ^|l,2;fc,A;), 
fc=i 



l,2;{A:,/}) = i{|l,2;fe,0 + |l,2;/,fe)}-^|l,2), 
l,2;[A;,/])=^{|l,2;fe,/)-|l,2;/,/c)}, 



d-2 



l,l,l;fc)=^|l,l,l;A;,/,0, 



(5.8) 

(5.9) 

(5.10) 

(5.11) 
(5.12) 

(5.13) 



1=1 
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\l,l,l;{k,l,j})=\l,l,l;k,l,j) 

-^|i,i,i;j)-%|i,i,i;0-%|i,i,i;fe). (5.14) 

add 

It follows from this classification of states that they will not mix with each other 
when the interactions are switched on, except perhaps for the two vector multiplets 
at n = 3 that are indistinguishable from the symmetry point of view. 

5.3 Eigenvalue equation 

For specified quantum number n and representation of 0(d — 2), the states listed in 
table 1 span a subspace of lowest-order energy eigenstates with definite symmetry 
properties. The corresponding first-order energy shifts E^ ■ are obtained from the 
matrix elements of the interaction 

S,nt = S2 + Ss (5.15) 

in these subspaces. 

Let 4'\{zq) and (j)^{zo) be the expressions in creation and annihilation operators 
at time zq that are associated to the basis vectors of a given subspace at level n. 
We label them by an arbitrary index A and define the normalization matrix Nab 
through 

{<PAmm)o = NABe-'''/^ (5.16) 

(we take t positive). The interaction matrix Has is similarly given by the correlation 
function 

{Mt)Sint<pUO))^ = {HAst + ■■■} e-"-*/^ (5.17) 

where the ellipsis stands for any terms independent of t. 

From the spectral representation of the correlation functions in the interacting 
theory, it is now straightforward to deduce that the solutions A of the generalized 
eigenvalue problem 

Hv = XNv (5.18) 

coincide with the first-order excitation energies 

AEl, = El,-El,. (5.19) 
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(a) 



(b) 



(c) 



(e) 



Fig. 1. Feynman diagrams contributing to the four-point function (5.20). In this 
figure time runs verticaUy from time 0, where the string modes are created, up to time 
t, where they are annihilated. 

If the representation of 0{d — 2) in the chosen subspace is irreducible, the eigenvalue 
problem is trivial, i.e. H is proportional to N and there is a single eigenvalue A. The 
n = 3 vector multiplets are an exceptional case from this point of view, where one 
has two eigenvalues that correspond to AE^ 2 ^^d ^-^3 3- 

5.4 Sample calculation 

For illustration let us consider the correlation function 



[a^.{t,p)ai {t,q)Sintaf{0,p)aj{0,q)) = {Mki,ijt -\ }e" 



-{p+q)t 



(5.20) 



By a decomposition into irreducible 0{d — 2) tensors, the interaction matrices Hab 
associated to the states in table 1 with two excited modes are then easily obtained 
from the matrix M^i^ij. Here we only discuss the computation of this matrix, which 
amounts to the evaluation of the 5 diagrams shown in fig. 1. 

The diagrams (a)-(d) are products of the free propagator (5.4) and the connected 
part of the correlation function 



(a-(t,p)5intar(0'P)>o = {^kit + •••}' 



-pt 



(5.21) 



Using Wick's rule the combinatorial factors are easily obtained from this expression. 
Moreover, from eqs. (5.2) and (5.3) one infers that a non-zero contribution to the 
diagram is obtained only when the interaction vertex is inserted at any time between 
and t. The part of the Feynman integral that adds to the matrix element M^i is 
then easily identified, and its value may finally be calculated following the lines of 
appendix A. In the case of the diagram (a), for example, this leads to the result 



M, 



(a) 
kl,ij 



Sk.iS, 



kiOlj ■ 



pq^ 



{d — 1) qr 



vr 



2C2 + 



{d+l)qr-^{d-l)'\cs] (5.22) 



at D = 2. 



16 



To evaluate the tree diagram (e), we only need to substitute the free propagators 
(5.2) and (5.3) for the lines from the sources to the interaction vertex. It again 
follows from these formulae that the non-zero contributions are obtained when the 
vertex is at times between and t. Some algebra then yields 

^kl,ij = P^^^'^ { i^ki^lj^Pq + ^kj^li + hl^ij (1 + 5pq)] C2 

+ [hi^lj (l + \5pq) + ^kjSli (5 + 6pq) + ^6klSij (1 + 6pq)] C3 | (5.23) 

for this diagram. 

5.5 Why is C3 = — 2c2 a special case? 

If the couplings are related in this way, there are some important simplifications in 
the Feynman rules that allow the matrix elements Hab to be calculated in closed 
form at any level n. We first observe that the interaction at D = 2 may be written 
in the form 

Sint = -\c2 I d^z {d+hd+h) (d-hd-h) , d± = do± idi, (5.24) 

if C3 = — 2c2. In tree diagrams the source-to- vertex propagators are thus 

{atizo,p)d±hi{w))^ = ±6kie{wo - zo)ipe-P^^"-'-^e^'P^\ (5.25) 

{a^{zo,p)d±hi{w))^ = ±5kieizo - u;o)^pe-P(^»--°)e^*P'"S (5.26) 

and using these expressions it can be shown that the diagrams vanish, except when 
the states connected to the vertex are as in the diagram (e) considered above. More- 
over, as can be seen from eq. (5.23), the contribution is then simply proportional to 
the identity matrix. 

The only other diagrams are tadpole diagrams of the type (a)-(d), which are also 
proportional to the identity matrix. Combining all contributions it turns out that 
the interaction is diagonalized by the Fock states |ni, . . . , Um', ki, . . . , km) and that 
the first-order energy shifts E"^ = [E^] are given by eq. (4.11). Whether there is a 
deeper reason for this (a hidden symmetry perhaps) is, however, not known to us. 
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6. Summary of results 

6. 1 Energy spectrum for general d, C2 and C3 

As explained above, the obvious symmetries of the effective theory do not exclude 
that the n = 3 vector multiplets mix with each other once the interactions are 
switched on. The detailed calculation shows, however, that this does not happen to 
the order of perturbation theory considered. For n < 3 the first-order energy shifts 
are then all of the form 

^h = Kl + 73 {^ [1^ (^ - 2) - ^] C2 + l^n,i (C3 + 2C2)} , (6.1) 

where the energy shift E^ ^ = [E^] of the ground state is given by eq. (4.9) and the 
coefficients i'n,i are listed in table 2f. 

Table 2. Values of the coefficients i'n,i in eq. (6.1) 



n,i 


^n,i 


n,i 


^n,i 


1,1 


^{5d + 7) 


3,2 


|(17d + 19) 


2,1 


^ (lid + 13) 


3,3 


i(9d + 43) 


2,2 


^ (lid + 13) 


3,4 


I {9d + 35) 


2,3 


^(5d + 25) 


3,5 


i(9d + 19) 


3,1 


i(17d + 19) 


3,6 


|(5d + 43) 



It follows from these results that most degeneracies of the energy levels are lifted 
by the interactions if C3 + 2c2 does not vanish. Equation (6.1) also shows that the 
first-order energy shifts tend to increase with n roughly like n^. The distances r at 
which the large-r expansion of the energy values applies must therefore be expected 
to scale with n. 

6.2 Energy spectrum in d = 3 dimensions 

Since h{z) has only one component in 3 dimensions, the interactions 5*2 and 5*3 are 
the same and all energy shifts E^ ■ are therefore proportional to C2 + C3. Without loss 



y One readily checks that the average energy shifts [-E^] obtained in section 4 from the partition 
function and the energy values quoted here satisfy the sum rule ^J 'Wn,iE^ ^ = Wn[E^]. 
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Fig. 2. String excitation energies /S.En = En — Eo in four dimensions [eq. (6.3)] for 
coupling C2 = (2(t)~ . Physical units are defined by the Sommer scale ro = 0.5 fm and 
we have assumed that TqCt = 1.36, as suggested by recent simulations of the SU(3) 
gauge theory. 

we may set C3 = — 2c2, for example, which shows at once that this is a case where 
the free-string degeneracies of the energy levels are not lifted, i.e. to this order of 
perturbation theory there is a single energy value E^ for each n. 

The duality constraint (3.9) now implies C2 = {2a)~^ and the expansion 



En 



TT 



ar + fi-\ [n- ITT 



TT 



2ot3 



n 



^) +0(r-^) 



(6.2) 



of the energy values to this order thus involves no unknown parameters other than 
the string tension a and the mass fi. Note that the "multiplets" listed in table 1 are 
either one-dimensional or do not occur in three dimensions. 

6.3 Energy spectrum in d = A dimensions 

Open-closed string duality requires C3 = — 2c2 in this case, so that there is again 
only a single energy value £"„ for each n. The coupling C2 remains unconstrained, 
however, and the energies 



Er, 



IT 



ar + fi-l — (« - ]^) 



C2 7j(" 



fe)' + 0( 



12 



(6.3) 



depend on this coupling in a non-trivial way. 

For illustration only, we plot the three lowest excitation energies in fig. 2, in units 
of vr/r and setting C2 = {2a)~^ . Taking this value as an example is suggested by 
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our results in three dimensions and also by the fact that the coupling derived from 
the classical Nambu-Goto action has this value (the relation C3 = — 2c2 holds in this 
theory too). The figure shows quite clearly that the corrections to the free-string 
energy splittings grow with n and that r must therefore be taken to increasingly 
large values to reach the asymptotic regime. 



7. Concluding remarks 

The calculations presented in this paper provide some insight into the effects of 
the string self-interactions on the energy spectrum in the effective string theory. In 
this approach the possible form of the interaction terms is constrained by symmetry 
considerations only, and the results are thus representative of a whole class of string 
theories rather than of any particular model. 

An important qualitative outcome in space-time dimensions d < 4 is that the 
free-string degeneracies of the energy levels are preserved to first order in the in- 
teractions. Since the number of interaction terms is rapidly increasing, it seems 
a bit unlikely, though, that the degeneracies will persist at higher orders of the 
perturbation expansion (unless the couplings assume some special values). 

Our calculations also show that the first-order correction to the ground-state en- 
ergy (the 1/r^ term) is significantly smaller than the corrections to the excited-state 
energies. Maybe this is why string behaviour in the static quark potential in SU(A^) 
gauge theories is seen to set in at surprisingly short distances [1], where the energy 
spectrum is otherwise still far from being string-like [9]. However, the question that 
will have to be answered first is whether the low-lying energy values in the gauge 
theory are indeed matched by the string-theory formulae (6.2) and (6.3) at large 
distances r. The issue can presumably be resolved in lattice gauge theory, using 
numerical simulations, although this will require a dedicated effort and maybe more 
powerful simulation techniques to be able to reach the relevant range of distances 
with controlled systematic errors. 

We are indebted to Hiroshi Suzuki and Wolfgang Lerche for their comments on 
the notion of open-closed string duality, and one of us (M.L.) wishes to thank the 
Yukawa Institute at Kyoto, where part of this paper was written, for the generous 
hospitality extended to him. 
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Appendix A. Dimensional regularization 

The effective string tlieory lives on a manifold with boundary where the standard 
momentum-space techniques of dimensional regularization do not apply. There is a 
more general formulation of dimensional regularization, however, that can be used 
here since it is based on the heat-kernel representation of the propagator in position 
space [26]. 

A.l Heat kernels 

The coordinates of a point z on the world-sheet are zq,zi, . . . , zo-i, where the last 
D — 2 parametrize the unphysical extra dimensions. We take these to be a product 
of circles, with circumference L, and now first need an explicit representation of the 
free massless propagator on this manifold, with Dirichlet boundary conditions at 
Zi = 0, r and periodic boundary conditions in all other directions. 

To this end we introduce two types of one-dimensional heat kernels, the first 

<(x,y;r) = ij;e^-"(-^)e-*-o, po = ^, m e Z, (A.l) 

Po 
for periodic boundary conditions and the second 

Kf {x,y;r) = - > sin(pia;) sin(piy)e~*^i, pi = , rii G N, (A. 2) 

Pi 

for Dirichlet boundary conditions. The heat kernel and the propagator on the world- 
sheet are then given by 

D-l 

Kt{z,w)=K^^{z^,w^-T)K^{zr,wr\T) J] Kf (z„,u;„; L), (A.3) 

a=2 
/■oo 

G{z,w) = / AtKt{z,w). (A.4) 

Evidently these formulae only make sense for non-negative integral values of Z) — 2, 
but an analytic continuation in the dimension of the world-sheet is possible later, 
after the positions of the vertices in the Feynman diagrams have been integrated 
over. 
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Using the Poisson summation formula, the one-dimensional heat kernels can also 
be written in the form 

oo 



n= — oo 



oo 
Kf(x,y;r) = (47rt)-l/2 ^ |g-(x-J/+2n.)V4t _ g-(a:+y+2nr)V4*| _ (^.6) 

n= — oo 

In particular, in the limit t ^ the sums 

Ko(t)=j;e-*Po, Ki(t) = j;e-*P^ (A.7) 

Pa Pi 

(where po and pi run over the values specified above) are equal to T{A7:t)~^''^ and 
— 1/2 + r(47rt)~"'^'^ respectively, up to exponentially small terms. 

A. 2 Calculation o/ (S'i)q 

Since the final results do not depend on the compactification scale L, we may now 
just as well set L = T in this subsection. When the heat-kernel representation (A. 4) 
of the propagator is inserted in eq. (2.13), the integral 

h f °° d 

{S,), = --{d-2)J^ dt [Ko{t)f-' -K,it) (A.8) 

is obtained, which is absolutely convergent for all complex D in the range HeD < 0. 
Moreover, by adding and subtracting the asymptotic term T(47rt)~^'^, raised to the 
power D — 1, from the first factor in the integrand, the expression can be analytically 
continued in D to the whole complex plane. This leads to the representation 

{Si)o = -^^ - 2) {l^ + /°° dt [^o(t) - r(47rt)-V2] ^^,(t)| . (A.9) 

at D = 2. 

We may now expand Ko{t) and Ki{t) in the exponentials e~^"^'^' '■** and e~y"-^/^) * 
respectively and integrate the terms of the double-series over the parameter t. The 
sum over m can then be performed analytically, and the formula 

^ n=l ■' 
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is thus obtained, which coincides with eq. (2.14). 

A. 3 More complicated diagrams 

The one-particle irreducible parts of the Feynman diagrams considered in this paper 
are tadpole diagrams where one or more loops are attached to a single vertex. All 
these diagrams can be worked out analytically along the lines explained above. 

In the first step the heat-kernel representation (A. 4) is substituted for all (internal) 
propagators and the integral over the world-sheet coordinates of the interaction 
vertex is performed. One is then left with an integral over the heat-kernel parameters 
that is absolutely convergent for sufficiently small ReD. The analytic continuation 
to the physical point D = 2 is finally achieved by separating the leading singularities 
of the integrand at small values of the heat-kernel parameters. 



Appendix B. Radon transform of spherically symmetric functions 

In this appendix we derive the expansion (3.2) from the spectral representation (3.1) 
and the fact that the Polyakov loop correlation function is spherically symmetric. 
The key point to note is that the left-hand side of eq. (3.1) has the form of an integral 
transformation that maps the correlation function to some function of z = \xi —yi\- 
This is actually a special case of the Radon transform [27], and we now first discuss 
some general properties of this transformation. 

B.l Definition 

In the following we consider smooth functions /(r) on the half-line r > that fall 
off faster than any inverse power of r at infinity. Their behaviour in the limit r ^ 
is left unspecified and can be arbitrarily singular. 

For any integer n > 2 and all z > 0, we define the Radon transform 9^„/(z) of a 
given function /(r) through 



fB.l) 



^nfiz)= / d"-iy/(r)| /3 



If we interpret r as the radius in M" , the integral on the right of this equation is taken 
over an ra — 1 dimensional hyperplane, parametrized by the coordinates y G M"-"^, 
at a distance z from the origin. It is easy to show that the transformed function is 

23 



again in the class of functions specified above, and the transformation can thus be 
apphed repeatedly. An immediate consequence of the definition (B.l) is then that 

In particular, 9^„ = 9^25^2 • • • 9^2 {n — 1 factors). 

B.2 Inversion of the Radon transform 
In n = 3 dimensions we have 

msf{z) = 27r drr/(r), (B.3) 

J z 

and the function /(r) can be reconstructed from its Radon transform through 

f{r) = -i27rr)-^-^9\3f{r). (B.4) 

The Radon transform 9^3 is thus an invertible mapping from the space of functions 
defined above onto itself. 

In view of the semi-group property (B.2), the same is also true for all other trans- 
formations. In particular, the inverse of 9^2 is obtained by first acting with 9^2 and 
then using eq. (B.4). Explicitly we have 

f (9^^^)2("-i) if n is odd, 
^-' = { (B.5) 

[(^3"^) ^"^2 if n is even, 

where ^3^f{r) = -{2TTr)-^f'{r). 

B.3 Example 

It is now almost trivial to find the function k{r) that satisfies 

yink{z) = e-^\ (B.6) 

If n is odd, say n = 2z/ + 1, it is given by 

— 

k{r) = {-2.r (^|:) e-^'^ = 2r (A) ' K..,{Er), (B.7) 
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where the normahzation of the Bessel functions is as in ref . [28] . In the other case 
the solution is obtained by noting that 

d\2f{z) = 2zK^{Ez) if /(r)=e-^^ (B.8) 

The apphcation of eq. (B.5) for n = 2u then yields 

— 

k{r) = {-2TT)-'' {^-^\ 2rK^{Er) = 2r(^y K^_^{Er), (B.9) 

i.e. the same expression in terms of the Bessel functions as above. 

B.4 Proof of eq. (3.2) 

As already mentioned at the beginning of this appendix, eq. (3.1) may be written 
in the form 

CXD 

9i,_i/(z) = j; |^;„|%-^"^ (B.IO) 

n=0 

where /(r) now stands for the Polyakov loop correlation function at distance r. The 
application of the inverse Radon transform and the results obtained in the previous 
subsection then immediately lead to the expansion (3.2). 

A technical point we would like to emphasize here is that the inverse of the Radon 
transform is unique and that this crucial property holds irrespective of the behaviour 
of the functions at small distances r. The Radon transform is actually local in the 
sense that the values of the transformed function at z > a only depend on the values 
of the input function at distances r > a, for any fixed a > 0. 



Appendix C. Open closed string duality and the Nambu Goto string 

The energy spectrum of a Nambu-Goto string with fixed ends, 

1/2 



f 27r 1 ' 

E^ = arll + —[-^id-2)+n]\ , n = 0,1,2, 



(C.l) 



was calculated some time ago by Arvis [29] through formal canonical quantization. 
Since the degeneracies of the energy eigenstates are the same as in the case of the 
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free string, the associated partition function reads 

oo 

Z = 5^u;„e-^"^. (C.2) 



n=0 



We now derive the closed-string representation of this partition function. 
Our starting point is the formula 

^=^[ W^ e-i^-(^o^)^* {g--r,(g)}'"', q = e'^^^^^*, (C.3) 

which is easily verified by expanding the curly bracket in powers of q and by inte- 
grating the resulting series term by term. The modular transformation property of 
the 77-function may now be used and this leads, after re-expanding the curly bracket, 
to the representation 

in which the closed-string energies are given by 

En = aT\l + -^[-±{d-2)+n]\ . (C.5) 



CjT2 

Integration of the terms in the sum (C.4) finally yields a series of Bessel functions 
of the form (3.2), with 

\Vn\ =---[-] Wn, (C.6) 

En ^^^ 

and open-closed string duality thus holds exactly. 



References 

[1] M. Liischer, P. Weisz, J. High Energy Phys. 07 (2002) 049 
[2] P. Majumdar, Nucl. Phys. B664 (2003) 213 

[3] M. Caselle, R. Fiore, F. Gliozzi, M. Hasenbusch, P. Provero, Nucl. Phys. B486 
(1997) 245 

26 



[4] M. Caselle, M. Panero, P. Provero, J. High Energy Phys. 0206 (2002) 061 

[5] M. Caselle, M. Hasenbusch, M. Panero, J. High Energy Phys. 0301 (2003) 057; 

Nucl. Phys. (Proc. Suppl.) 129-130 (2004) 593 
[6] M. Caselle, M. P epe, A. Rago, Nucl. Phys. (Proc. Suppl.) 129-130 (2004) 721; 

|hep-lat/0406008| 



[7] M. Caselle, M. Hasenbusch, M. Panero, |hep-lat/0403004 

[8] S. Perantonis, C. Michael, Nucl. Phys. B347 (1990) 854 

[9] K. J. Juge, J. Kuti, C. J. Morningstar, Nucl. Phys. B (Proc.Suppl.) 63 (1998) 

326; ibid. 73 (1999) 590; ibid. 83 (2000) 503; ibid. 106 (2002) 691; Phys. Rev. 

Lett. 90 (2003) 161601 
[10] G. Parisi, R. Petronzio, F. Rapuano, Phys. Lett. B128 (1983) 418 
[11] P. de Forcrand, G. Schierholz, H. Schneider, M. Teper, Phys. Lett. B160 (1985) 

137 
[12] S. Ohta, M. Fukugita, A. Ukawa, Phys. Lett. B173 (1986) 15 
[13] B. Lucini, M. Teper, Phys. Rev. D64 (2001) 105019 
[14] K. J. Juge, J. Kuti, F. Maresca, C. Morningstar, M. J. Peardon, Nucl. Phys. 

(Proc. Suppl.) 129-130 (2004) 703 
[15] Y. Nambu, Phys. Lett. B80 (1979) 372 

[16] M. Liischer, K. Symanzik, P. Weisz, Nucl. Phys. B173 (1980) 365 
[17] M. Liischer, Nucl. Phys. B180 (1981) 317 
[18] K. Symanzik, Nucl. Phys. B190 [FS3] (1981) 1 
[19] M. Liischer, Nucl. Phys. B254 (1985) 52 
[20] K. Dietz, T. Filk, Phys. Rev. D27 (1983) 2944 
[21] J. Ambj0rn, P. Olesen, C. Peterson, Phys. Lett. B142 (1984) 410; Nucl. Phys. 

B244 (1984) 262 
[22] M. Liischer, Commun. Math. Phys. 54 (1977) 283 
[23] K. Osterwalder, E. Seller, Ann. Phys. (NY) 110 (1978) 440 
[24] E. Seller, Gauge theories as a problem of constructive quantum field theory and 

statistical mechanics. Lecture Notes in Physics 159 (Springer, Berlin, 1982) 
[25] M. Liischer, Selected topics in lattice field theory. Lectures given at Les Houches 

(1988), in: Fields, strings and critical phenomena, eds. E. Brezin, J. Zinn-Justin 

(North-Holland, Amsterdam, 1989) 
[26] M. Liischer, Ann. Phys. 142 (1982) 359 

[27] I. M. Gelfand, M. I. Graev, N. Y. Vilenkin, Generalized functions. Vol. 5 (Aca- 
demic Press, New York, 1966) 
[28] I. S. Gradshteyn, I. M. Ryzhik, Table of integrals, series and products (Academic 

Press, New York, 1965) 
[29] J. F. Arvis, Phys. Lett. B127 (1983) 106 

27 



